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Abstract 

It is shown that gauge theories with fermions are most naturally stud- 
ied via a polar decomposition of the field variable. This is the fermionic 
analog of the preprint cond-mat/02 10673. The hope is that these two 
put together will enable the treatment of neutral nonrelativisitc matter 
composed of electrons and nuclei in a nonperturbative manner with nuclei 
and electrons treated on an equal footing. We recast the electron-phonon 
(superconductivity) problem in the hydrodynamic language and indicate 
how it is solved. In particular we focus on the a.c. conductivity. 

1 Introduction 

Read My Lips ! 

George Bush Sr., 1993 Campaign. 

The program of quantizing hydrodynamics has a long and distinguished his- 
tory. Landau [1] and his students were among the first to attempt this. Later on 
Sunakawa et.al. [2] and others - notably Rajagopal and Crest [3] took this pro- 
gram further. Dashen, Sharp, Menikoff and Goldin[4] in the seventies introduced 
many of these ideas. Recently, Jackiw and collaborators [5] have revived interest 
in this approach in the context of relativistic quarks. In our earlier work, we 
introduced the DPVA for fermions[6]. We also note that Rajagopal and Grest[3] 
had already in the seventies pointed out the need for having a nonzero-phase 
functional found in the DPVA. In our earlier work[6] we made a first pass at 
computing the phase functional. This attempt yielded an answer that in retro- 
spect is quite wrong. Upon closer examination the Uo{q) of our earlier work[6] 
is imaginary when it was postulated to be real(for small q). So far the au- 
thor has avoided this issue by taking refuge under the the sea-boson approach 
that enables us to derive the momentum distribution, anomalous exponents, 
quasiparticle residue and so on without yielding the full dynamical propagator 
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which is of interest only because it contains information about quantities just 
mentioned. If one is able to compute them without having to compute the full 
propagator so much the better. However, there are physical problems in which 
the full propagator is important. The X-ray edge problem [8] is one such. In fact 
wc tried using the DPVA to compute the X-ray edge spectra in a preprint [7] 
and found that we obtain the right answers in one dimension but the answers 
in higher dimensions were inconsistent with Mahan's exact results [8]. Thus we 
must now face up to the harsh reality and try and address this (hopefully last) 
hurdle. This article is the fermionic analog of the preprint cond-mat/0210673. 
In an earlier preprint, after much reflection, we chose to dismiss the approach 
that only uses the hydrodynamic variables namely the density and its conjugate 
as 'myopic' (mypoic bosonization). This is beacause a hamiltonian formula- 
tion in terms of the hydrodynamic variables is unable to distinguish between 
fcrmions and bosons. We have to further decompose these variables in terms of 
linear combination of oscillators in order to distinguish between the two statis- 
tics. However, the sea-boson approach is not without its share of problems. 
For one it does not generalise to finite temperatures easily. Also the full dy- 
namical propagator is not reducible to quadratures due to a technical difficulty. 
Both these problems may be resolved in an approach that incorporates only 
the hydrodynamical variables. We show in this preprint, that the path integral 
approach is an avenue to distinguish between the statistics when using only the 
hydrodynamical variables. 

The author had this idea in 1993 but conversations with various knowledge- 
able people suggested that this approach that only uses hydrodynamic variables 
is not going to work out for fermions, since one has to take into account the 
extended nature of the Fermi surface. However, this idea seems too important 
to pass up. In particular, the natural manner in which gauge theory may be 
studied in this approach makes this effort for fermions worthwhile and urgent. 

2 The Field Operator 

In our earlier work [6], we showed how the field operator may be expressed in 
terms of currents and densities. We reproduce the formula here. 

V>(x) = e - in ( x >e < *(M ;x >Vp(30 (1) 

The conjugate n obeys canonical commutation rules. 

[H(x),H(x)] = (2) 

[n(x),p(x)]=z«5 d (x-x) (3) 
[p(x),p(x')]=0 (4) 

There are some technical difficulties associated with the fact that n is not strictly 
self-adjoint, but we shall operate in the high enough density limit where we may 
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assume that we expand around a nonzero mean for the density and this means 
that we may choose II to be self-adjoint. The conjugate II may be related to 
the current as follows. 

J(x) = -p(x)Vn(x) +/9 (x)V$([p];x) -p(x)H<&(H;x),VII(x)] (5) 

The phase functional $ determines the statistics of the field tp. As shown 
earlier [6] $ obeys a recursion relation as depicted below. It must be pointed out 
that the need for having a nonzero $ for fermions was appreciated by Rajagopal 
and Grest[3] way back in the seventies. However, the constraint below brought 
about by the imposition of Fermi statistics on t/j was probably first shown in 
our earlier work [6]. 

%(y)-^(y-x)}];x)-$(H;x) 

-$([{ /9 (y)-,5 d (y-x)}];x) + cI.([ /5 ];x')-7rm(x,x') (6) 

for x 7^ x' and m(x, x') = — m(x',x) is an odd integer. A further constraint 
on <E> emerges when we realise that Eq.( 5) has to be consistent with current 
algebra. We know that currents and densities obey the current algebra. In other 
words, 

[J(x),p(x)] =-zp(x)V x <5 d (x-x) (7) 
[J 4 (x),J,(x')] = -iJ J -(x)Vj c <J d (x-x' ) +U l (x')V J x ,(5 d (x' -x) (8) 

[p(x),p(x')]=0 (9) 

In what follows we try and impress upon the reader the need for a systematic and 
general approach by pointing out that reasonable sounding and rather general 
ansatzs fail to obey the constraints just outlined. Thus it seems that constraint 
on $ due to current algebra is almost inconsistent with the constraint due to 
the statistics requirement. 

2.1 A Serendipitous Surmise 

Many years ago the author had a conversation with then the student now Prof. 
A.H. Castro-Neto where the latter suggested that maybe the field operator is 
simply given by, 

V>(x) w e _i ^ e,q XXq Vp° (10) 

where by definition Xq = iq • j(— q)/(q 2 A^°) and p° = N°/V. Later he realised 
that if we choose LY q , X i] — as is in fact mandatory 1 , then fermion commu- 
tation rules are not obeyed. However one may take solace in the fact that at 
least one commutation rule does come out right namely [ip(x), p q ] = e* q x ^(x). 

1 Since strictly speaking it is the conjugate to p as denned by the line integral of the ratio 
of the current and density that enters; these commute amongst themselves. 
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A refinement over this ansatz was attempted in our earlier work[6] by introduc- 
ing an additional phase functional of the density linear in the density and this 
was also inadequate since by now the author knows that $ there when com- 
puted was imaginary when it was postulated to be real. A compromise was also 
suggested that involved multiplying and dividing by the free propagator and 
using the exact version in the numerator and the bosonized free propagator in 
the denominator. This trick though repugnant to most, gives us an anomalous 
exponent of the Luttinger liquid as we shall see below. However, this anomalous 
exponent is off by a factor of two from the exact one obtained by Mattis and 
Lieb[9]. From Eq.( 10) we may write, 

G(x-x ) = ^ f (x')V(x)J) « ^ e ^ q ( eq ~ E<1 ) Xq ^ p° 
Again using the trick outlined in our earlier work[6] we may write, 




Here Go(x — x ) is the propagator obtained from elementary considerations. In 
one dimension, we may see that (A q A_ q ) w k F S(q)/(q 2 N°). The structure 
factor S (q) = \q\/(2k F ) for the interacting case. For the interacting case we 
have, S(q) = (v F /v eff ) S (q). 

G(x-x ) = G Q (x-x ) e Jo M \ 2v '" 2 J 

~Gb(*-z') (^) 7 (13) 

where 7 = ^TJJ ~~ \ ■ This exponent is exactly one half of the exponent obtained 
by Mattis and Lieb[9]. What is even worse is, we have shown in an earlier 
preprint [7] that when applied to the X-ray edge problem, we obtain the well- 
known results of Mahan in one dimension but not in higher dimensions. Thus 
it would appear that there is something amiss in the expression for the field 
operator. The present attempt is to finally address this difficulty. 

2.2 No-Go Theorems 

We make some observations that render seemingly obvious and promising ap- 
proaches futile. 

The First No-Go Theorem : 
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Here we show that a simple and very reasonable ansatz for m(x, x ) in one 
dimension, fails. Set m(x,x ) = 9(x — x ) — 6{x — x). Clearly m = ±1 is odd 
and m{x,x ) = —m(x 1 x). It may be shown using the method of generating 
functions that the most general solution to Eq.( 6) is as shown below. 

*(\p];x)=w [ X dy p(y) + C ([{p(y)-5(y-x)}})-C ([p]) (14) 

J — oc 

Here Co is arbitrary. Unfortunately, the presence of the first term means that 
current algebra is violated. Using Eq.( 14) in Eq.( 5) we have 

J(x) = -p(x)d x U(x) + p(x)np(x-) (15) 

Then the current-current commutator reads as follows. 

[J(x), J(x )] = — iJ(x)d x 5(x — x ) + iJ(x )d x >5(x — x) 

—inp(x)p(x-) d x S(x — x ) + inp(x )p{x_) d x 'S(x — x) (16) 

The last two terms tell us that something is not right. Thus this attempt has 
come to naught. It is telling us that perhaps m depends on the density as well. 

The Second No-Go Theorem : 

Here we show that a general looking ansatz for $ that manifestly obeys current 
algebra fails to reproduce fermion commutation rules. 

$([p];x) = B ([{p(y) 5 d (y x)}]) - Bl ([p]) (17) 

Using the fact that V B ([{p(y)-S d (y- X )}}) = {-iB ([{ P (y)-5 d (y-x)}}), Vn(x)] 
(this can be shown easily by Fourier decomposing <J> with respect to p). 

J(x) = -p(x)Vfi(x) (18) 

n(x) = n(x) + [iB 1 ([ P ]),n(x)] (19) 

It can be seen that II is also a canonical conjugate of p. 

[H(x),H(x)]=0 (20) 

[n(x),p(x')]=^ d (x-x') (21) 

Thus Eq.( 18) obeys current algebra. Unfortunately Eq.( 17) fails to obey 
Eq.( 6). To see this we merely plug in Eq.( 17) into Eq.( 6) and find, 

F([p];x) - F([p];x') =? tt (odd integer) (22) 

F([p];x) = B.dipiy) S d (y x)}]) - B ([{p(y) 5 d (y x)}]) (23) 
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It can be seen that Eq.( 22) is a contradiciton. If F([p];x)/n is odd then 
F([p];x )/tt has to be even for all x ^ x. In particular, F ([p\; xl^ / n and 
x 2 )/7r have to be both even. But now their difference can no longer be odd, 
thus violating the recursion relation. While this is somewhat disappointing, it 
is not surprising since Eq.( 18) suggests that the velocity is irrotational, namely, 
v = —VII. This is clearly not general enough. We would like to be able to write 
v = —VII + v r . In what follows we shall do precisely this. 



We set J ~ p v and ask what properties should v possess ? This implies the 
following commutation rules on the velocity operator. 

[v(x),p(x')] = -iV x <5 d (x-x') (24) 

and, 

-ip(x)6 d (x - x'jViv^x) + ip(x')5 d (x - x )V J x ,v,(x ) 

+p(x)p(x')[v i (x),v J (x')] = (25) 

From Eq.( 5) it is clear that we may expect to be able to write the velocity as a 
sum of two parts, one irrotational and the other depending only on the density 
which has a nonvanishing curl. Thus, 

v = v s + v r (26) 

Here v s = —VII s for some II S . Without loss of generality we may choose 
[v s (x), v|(x )] = 0. In other words we include only those parts of the velocity 
into the irrotational part that makes it commute amongst itself. Of course, 
[v[(x),vj(x )] = 0. Crucially however, [v[(x),v|(x )] ^ 0. Indeed, 

K(x),v s (x )] = -V^K(x),^] (27) 

[v s (x),p(x)] = -zV x <5 d (x-x') (28) 

and, 

-zp(x)S d (x - x')VivJ(x) + ip(x)5 d (x - x )V^v[(x ) 

-p(x)p(x')Vi[i^,v;(x')]-p(x)p(x')Vi4vr(x),i^]=0 (29) 
This means that, 

~ ^' ( p(x) X } V ^ (X) ~ V ^^x)'^ (X )] = fM>**' ) ( 3 °) 

for some /y([p];x;x ) = /^([pjjx ;x). Therefore so long as Eq.( 30) is obeyed 
then current algebra is also respected. In light of the two no-go theorems let us 
set $ to be as shown below. 

$(M;x) = * <rr (W;x) + B ([{p(y) - S d (y - x)}]) - B ([p\) (31) 
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Here <&irr stands for some particular solution and <I> would be a general solution. 
If $i rr obeys current algebra and the recursion relation so too does From 
Eq.( 5) we may read off the following formula for the velocity operator v = 
v s + v r . Where, 

v'Xx) = V$ irr (W;x) - [$ irr (W;x), V x j^—] (32) 

and v s = — VII s where II s is given by Eq.( 19). We are looking for a singular 
solution to $ that also obeys the statistics requirement. Anticipating some 
simplifications, let us reparametrise <& as follows. 

<fr irr ([p];x) = vr L([{p(y) -S d (y- x)}];x) - tt L([p];x) (33) 

In fact it can be shown that only a form such as in Eq.( 33) is consistent with 
the statistics requirement in Eq.( 6). The statistics requirement may now be 
written in exponential form as follows. 

e - 4 7rL([{p(y)-,5' i (y-x')}];x) e47 rL([p];x) el 7 r L([{p(y)-^(y-x)}];x') e - 47 rL([p];x') 
= _ e - 4 7rL([{p(y)-,5' i (y-x')-5 ti (y-x)}];x) et7r L([{p(y)-^(y-x)}];x) 
x ^L([{p(y)-5 d (y-x')-,5 d (y-x)}];x') e - 4 7rL([{p(y)- I 5 £i (y-x')}];x') (34) 

The requirement in Eq.( 34) is unfortunately strictly speaking, inconsistent with 
the constraint from current algebra. To see this, we note that in one dimension, 
the constraint from current algebra reads as follows. 

'- ! ^' /(l)HH ^' /W1 ^ 
This can be satisfied only if, 

v r {x) = -ir[d XJ ^r,L(\p])] (36) 
dp(x) 

for some L (see appendix A). We see that this means that the L of Eq.( 33) has 
to be now independent of x. This means unfortunately that Eq.( 34) is violated. 
However, one can take the point of view that we would like the current algebra 
to be 'almost obeyed', in other words we allow L to depend extremely weakly 
on x. This will also open up the possibility that Eq.( 34) may now be satisfied. 



3 Evaluating the Phase Functional 

Here we would like to compute the phase functional that obeys the constraints 
just described. Computing the phase functional systematically seems particu- 
larly difficult. Intuition suggests that the phase functional ought to be related 
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to the properties of the free theory. Specifically, we should be able expand it 
in powers of the density fluctuations and somehow relate the coefficients to the 
density-density, current-density and current-current correlation functions. From 
the preceeding sections, it appears that we may write, 

V>(x) = e lA([p] ' x) e- m( ^^(xj (37) 

and A is some functional that depends weakly on x. In fact it has to be strictly 
independent of x in order for current algebra to be obeyed. However, the weak 
dependence is needed to recover Fermi statistics. Indeed if we choose A to be 
independent of x and then evaluate the equal-time version of the propagator of 
the Luttingcr model, we find an anomalous exponent equal to one-half of the 
exact one derived by Mattis and Lieb. In other words, an A independent of x is 
more or less equivalent to the serendipitous surmise. Thus we have to do better. 
In general II may be written as IT(x) = X + fl(x). Here X is conjugate to 
the total number of particles and II is strictly selfadjoint since it is related to 
currents and densities. We shall assume that X is also self-adjoint but may 
not be expressed in terms of Fermi bilinears. It may be shown[6] that X has 
an expression shown below. 

X ° = E np M Ln{ c k ) - ^ Y, n ^( k ) L < c l ) ( 38 ) 

k k 

It seems that the lagrangian formulation is appropriate here. The hamiltonian 
formulation has not worked out for unknown reasons. In the lagrangian formu- 
lation, all the variables are c-numbers. In the usual path integral formulation, 
the Fermi fields are complex Grassmann numbers and not ordinary complex 
numbers. As pointed out earlier, we take the point of view that a complex 
Grassmann number field may be 'simulated' by a complex ordinary number 
field provided one introduces an appropriate phase functional A. The value of 
A is tuned so as to reproduce the correlation functions of the free theory. The 
anticommutating nature of the Grassmann number is not present in the polar 
decomposition in Eq.( 37) but is recovered at the level of the propagator. To 
see this we note that in the path integral formulation the fields obey periodic 
(KMS) boundary conditions [11] namely, 

V>(x,i-i/3) = -e^ V(x,t) (39) 

Using Eq.( 39) in Eq.( 38) we find, 

X (t - i/3) = X (t) + i Log (-e^) = X Q (t) - tt + i/3/i (40) 

That is, the one-particle propagators obey Fermi commutation rules mainly due 
to the fermionic boundary conditions obeyed by the global Klein factor. The 
current operator may be written as follows. 

J(x) = -p(x)Vn(x)+p(x)V([{p q - e ^ x }];x) (41) 
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where, 

VA([p];x)=V([p];x) (42) 

First we note that in order for current algebra to be obeyed, A ([/?]; x) has to be 
independent of x. In other words V([p];x) = 0. Unfortunately this means that 
Fermi statistics is no longer obeyed. However, we can take the point of view 
that A depends very weakly on x and thus we may ignore the derivative of A 
since it involves taking the difference between neighboring points. On the other 
hand, for statistics we are usually interested in the opposite extreme namely 
the asymptotic regime |x — x | — > oo. In other words, we may no longer ignore 
A([p];x) — A([p];x ) for large separations. However ignoring V also seems to 
be a bad idea. For example if we do and use the Dashen-Sharp formula for the 
kinetic energy we get a hamiltonian that describes bosons rather than fermions. 



K 



' d d x ( J 2 (VpY 



2m \ p 



Thus it appears that we have to retain V as well. We would like to write down 
the action in the lagrangian formulation. For this we first write, 

I= /*U.V,. t ^ (44) 



2m V 4/? 
where, 

v'([p];x)=V([{p q - e ^ x }];x) (45) 

Till now the discussion has been at the quantum level. Now we would like to 
recast this theory in the path integral language. For this we have to compute the 
classical action in terms of the collective coordinates n, p. This means we have 
to 'de-quantize' the quantum hamiltonian in Eq.( 44). Happily, while there are 
many inequivalent ways of quantizing a given classical hamiltonian there is likely 
a unique way of de-quantizing a quantum hamiltonian. Thus we may make the 
following observations. The Hamilton equations for the canonical variables are, 

d t p( X ,t) = -^-V x (2p(x)(V x n(x) - V ([/>]; x))) (48) 
L(p, d t p) = J d d x n(x, t)$p(x, t) - H(p, n) (49) 
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L( P ,d tP ) = -J d d xn( x ,t)^-v x (2 P (x)(v x n(x)-v'(W;x)))-y ^^(vn-v 

= ^Jd d x V x n(x,t) ( /5 (x)(V x n(x)-V'([p];x)))-|^^(Vn-V') 2 + ^ 



= ^Jd d x p(x) V'([p];x) (V x n(x)-V'([p];x))+ 1 ^ (p(VIT- V) 2 - ( ^ '" 



4p 
(50) 



4 One Dimension 

In one dimension we have, 

d t p(x,t) = ~d x (p(x)(d x IL(x) - v'([p};x))) (51) 

(d x U(x,t)-V'([p];x))=-^— f dxd t p(x,t) (52) 
The Lagrangian is given by, 

/OC r-X 
dx V'([p t ];x) / dx' d t p(x ,t) 
-OC J —oo 

f°° dx ( m 2 / r , , , , A 2 (d x p(a;,*)) 2 \ 
The current operator is then given simply by, 

r 

J(x, t) = ml dx d t p(x , t) (54) 

J — oo 

So lar the discussion has been exact. Now we would like to make some approx- 
imations. In particular we assume that we are in the degenerate regime so that 
p(x,t) = po + p(x,t) and p(x,t) <C po and p — N a /L the number of electrons 
per unit length. 

/OO f-X 
dx V'([p +pt];x) / dx d t p(x ,t) 
-oo J — oo 

f°° dx (m 2 ( { x , , A 2 (d x p(x,t)) 2 \ 
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P(a,t) = |5>,(t)e-** (56) 

/oo /-a; 
dxV'([p +p t ];x) / dx d t p(x ,t) 
-oo J — oo 

/oo j 2 / ^a; \ 2 



9#0 



dx <9 t p(x , t) = j ^ dtP q (t) — 
Hp, dtp) = ^v'([p + pt];q) 

7 777 

j,(t) = — (59) 



«'([/>o + &];?) = f ^V(lpo+Pthx) e-** 



(60) 



We would now like to determine w and hence A by forcing this action to repro- 
duce the correct current-current, current-density and density-density correlation 
functions of the nonintcracting Fermi theory. Unfortunately since the time de- 
pendence of v is through p alone this is not possible as we have found after 
repeated attempts by choosing various v , linear, quadratic and so on in p. Thus 
we now have to form another reinterpretation. We have the luxury of reinter- 
preting these formulas anyway we wish so long as the properties of the free 
theory are properly recovered. Thus we shall transfer the time dependence on 
to the function A itself thereby making the action depend on the history of the 
density of the system. 



5 Lagrangian Formulation 

The discussion in the above sections shows that for fermions, starting from a 
Hamiltonian formulation and then moving to the Lagrangian formulation does 
not work out. Thus we would like to start with a Lagrangian formulation at the 
outset. For this we postulate that the field operator is given by, 

V>(xi) = e * A (M; x *) e-^^VpR) (61) 
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The notation A([p];xi) is meant to imply that the phase functional potentially 
depends on the history of the density configurations of the system. It does not 
imply that the phase functional is explicitly time dependent. Also A depends 
extremely weakly on x thus we may ignore VA w 0. However as pointed out 
before, we may not ignore A([p];xt) — A([p];x t) for |x — x | — ► oo. All the 
variables are c-nmbcrs in the Lagrangian formulation. Even the field variable 
in Eq.( 61) is a c- number. The anticommuting nature of the field variable is 
captured in the global Klein factor as already pointed out earlier. The current 
is now given in terms of the canonical variables as follows. Since we may ignore 
the gradient of A we have 

J(xt) = -p(xi)Vn(xt) (62) 
The action may now be recast as follows. 

d d x 



Sfree = J dt J 



pd t U - pd t A - p(WU) 2 {VpY 



if3 

" '^J (63 » 

After this reinterpretation we would now like the action to be purely quadratic 
in the canonical variables and yet be able to reproduce the correct fermion 
correlation functions. Thus we make the following assumption. In the high 
density limit we may ignore the density fluctuations whenever appropriate and 
this is equivalent to working in the asymptotic regime. 

p(xt),^ £ -^ qlie -^ (64) 



qn 



II(xt) = ^e^ x AVe z "* (65) 

qn 

A([p];xi)= J2 e iq - x A([p];qn) e^* (66) 

q#0n 

The Klein factor which is the position in dependent part of II ensures that the 
fermion KMS boundary conditions are obeyed by the field variable since we still 
have z n — 2irn/ (3. Here p is a real function hence [p(x, t)\* = p(x, t*) = p(x, —t) 
since t € [0, — i0\. This means p* „ = p-q, n - For a similar reason A* n = A_ q .„. 
Note however that n is not self adjoint due to the first (q = 0) term, namely the 
non-self adjoint Klein factor. The position dependent part of n is self adjoint. 
Also we have [A([p]; q, n)]* = X([p\; — q, n). 

The action in Eq.( 63) would be identical to the bosonic action were it 
not for the term — J pdtA. Thus this term is crucial and leads to the fermionic 
correlation functions for the current-current, current-density and density-density 
correlations since the current is still given simply as in the bosonic case. The 
action may now be written similar to the bosonic case, 

Sfree = ^(-i/3Zn)Pqn^qn + i(3N° ^ q 2 A q „A_ q _„ 
qn qn 
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qn q/On 

where p° = N°/V is the mean density of fermions. We have included the leading 
anharmonic corrections since this is going to be important in future works where 
we may choose to go beyond RPA. We choose, 

H[p]i W) = C(qn) p-q- n + S nfi D([p];q) (68) 

We would like to compute the unknowns C and D in the above ansatz by forcing 
this action to recover the right dynamical current-current density-density and 
current-density correlation functions of the free Fermi theory. 



Sfree = ^H/3z n ) p q „X q „ + i(3N° ^ ^X^X. 



q,— n 



+ TK[0 Y qVq,np-q,-n + if3 ^ Z n C(qn) PqnP-q ,_„ (69) 



4N a ■ 

qn q^On 



We may now integrate out X to arrive at an effective action, 



qn 



+ 4^0 Y ^Pq.nP-q + ^ Z ™ ^fa") PwP-q-n (70) 

qn q^On 

Thus we have, 

From this, we may read off a formula for C in terms of the density-density 
correlation function we already know how to compute. 

In one dimension we may write, 

where vf is the Fermi velocity. Since we taken care to ensure that current 
algebra is obeyed we may expect to find that the current-current, current-density 
and current-current correlation functions are also properly recovered. The above 
action does not give everything right. For example, the three body correlations, 
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^Pq,nP_q+ q ' - n +n P-q -n J arc zcr0 f rom the above action but clearly they 
are nonzero in the full free Fermi theory. This means we have to go beyond 
the quadratic action to capture three and higher body correlations. This is not 
important in the RPA sense and shall ignore this. 

Finally, we have to make sure we are able to compute the full propagator. 
Unfortunately a straightforward application of the formula for A does not yield 
the right propagator. Thus we have to redefine the field as follows. We may 
decompose the field variable into fast and slow modes. It is the slow modes 
that contribute to the action. But the fast modes are needed in the one-particle 
propagator. Thus we write, ip = ^ fast ipslow We postulate that the fast modes 
are not affected by interactions, since the action only involves hydrodynamic 
degrees of freedom, in other words the slow modes. In other words, 



V>(xi) = */ ost (xi) e 



e iq.x e ,„t (Xqn _ G(qn) fi _ 







(74) 



Here G(qn) has to be recomputed since we are unable to make contact with 
the free propagator otherwise. Perhaps this suggests that the real theory will 
involve some very complicated forms of A. The main aim of this preprint is to 
convince the reader of the urgency of somewhow making this scheme work for 
then wc will be able to compute the observable properties of a theory of neutral 
matter where the only adjustable parameters are the relative concentrations 
of the various elements that make up the material. Now we specialise to one 
dimension where the analysis is somewhat simpler. In furture versions of this 
preprint, we address the 3d case and apply it to the electron-phonon problem 
briefly discussed in the next section. We would like to write the propagator in a 
form that may be decomposed into a product of fast and slow propagators. To 
this end we demand that only the resonant regime of the propagator be properly 
recovered. In other words when (x — x ) ~ ±up(f — t ). 



i)Ux't) Mxt)) = — e 4 ^*-*') 
/ 2ni 



D ikF [{x — X )—VF{t — t )] 



-ikp [(x—x )-\-VF{t—t )] 



(x — x) — Vp{t — t') (x — x) + VF(t — t') 



= ( 2i ) e -^Ut-t) 



sin[kp{x — x )] 
{x-x' f -v 2 F {t-t') 2 



In other words, 

** fast (x't) * f ast(xt)) = p° e - ife ?"(*-*') sin[k F \x-x'\] 
Thus we may deduce, 

2 



\x — X 



k F 



(x-x') 2 - v%(t-t') 



(75) 



(76) 



(77) 



14 



In what follows we shall only insist on being able get the resonant part of the 
slow propagator correctly. In other words, [(x — x ) 2 — v F (t — t ) 2 ] -1 . We shall 
not be too concerned if we fail to get the \x — x | in the numerator. From the 
hydrodynamic approach we have, 



2 /-*i q ^o, n y 



q (x-x ).Z„(t-l )_„»<! (X -•),»„( 



(78) 



F (q, n ) = (( x q,n - G(q, n) p- q -n) {X- q - n - G(-q, -n) p q . n )) 

= {X q!n X_ q _ n )-G(q,n) (p_ q _ n X_ q _ n )-G(-q,-n) (X q , nPq , n )+G(q,n) G(-q,-n) (p q . n p- q ,-n) 

(79) 

In one dimension we may deduce, 



If we set, 



then, 



(Pq,nP-q,-n) 



2N" z 







{■Xq,n-X—q,—n) 



(X q , n P q ,n) — ( a ■ 



2/3N \vj 

z 



z. 

Jq 



Zn +q 2 
F 

2 
n_ 

2 



F(q,n) 



2v 2 F 
f3N° 



(80) 

(81) 
(82) 

(83) 
(84) 



Here we have taken the liberty to ignore the pole at z n — since this may be can- 
celled by an appropriate time-independent additive contribution to A([p];x, t). 
This does not affect the action but serves to regularise the above integrals. 



Vslow( X 1 )^slow{xt) 



q (x-x ).Z„(t-l )_„iq (x -x) i n (t 



F(q.n) 

(85) 
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- (2-e 19 ( x - x ') e z " ( *-*' ) -e l<? {*'-*) e *n(t'-t)^ F(q,n) 

J dq ^iq (x-x') e i v F \q\{t-t) _ & iq (x -x) e i v F \q\{t-t')^ ^gg-j 



Thus we have, 



(87) 



Now we compute the total momentum-total momentum correlation function of 
the free Fermi theory. This will be important when we compute a.c. conductiv- 
ity. To this end we write, 



P(i) = - E (iq)p q , n X q y e (V 



■z„) t 



(88) 



q,n,n 



In one dimension, in the hydrodynamic language we may write, 
<TP(t)P(t')>= £ E {Pi,nX qtn .p q . tm X q . tn .) e K„' -*->*>. 



= 2 g 2 (p 9 ,„p_ g ,_„)(x ?y X_ gi _„,) e (V "*»)(*-*) 

- £ <z 2 (*, y p,y) c (v—) («-*') 

= E « 4 ^(^ + 4? 2 )" 1 ((^ + ^) 2 +45 2 )" 1 (*"*') 



q,n,m ^ ' 



EM 



/? 2 

1 



/J zl+4v 2 F q 2 



z m (t-t ) 



In the Fermi language, 



(3 zl + 4v 2 F q* 



„Z m (t-t ) = Q 

(89) 



kp 



T P(t) P(t)) = J2(k.k) (n k n k ,) = ]T fc 2 l fc D =N°f (90) 

k,k' k 

Thus there seems to be a discrepency. We must then take the point of view 
that the hydrodynamic approach predicts the right time-dependent part of the 
total-momentum total momentum correlation function. 



16 



6 The Luttinger Model 

The action for the Luttinger model reads as follows. 

Shutt = ^2(-i(3z n )p qin X qiTl + i[3N° q 2 X 5i „X. 5 . n 



q,n 



q,n 



-^o « P^P- q ,- n + *^ 2^ ijvo Pi> n P-i>- n + 



q^On 

(91) 



Here V^) > is a positive constant signifying repulsion. Thus we may write 

y eff — v f + 2 - Unfortunately it seems that we have to redefine G 
(we may euphemistically call this 'renormalization'). Thus we retry G(q,n) = 



\ — 



N° 



F(q,n) 



V 2 

V eff 



2?J 2 

X- 



2/37V /3N° 



AX 



2/3N° 



12, 2 2\~ 1 
(*n + v eff1 ) 



(92) 



I ^2 ^2 - e iq ( x ~ x ' ) e z n(t-t ) — e *9 (x'-x) e z„(t'-t)^ 



Jo q K 



£ iq {x-x') e iv eSS \q\(t-t) _ £ iq (x -x) e iv ef s \q\(t-t 



V -fl-X^- + . 

Av F V e ff 

(93) 



It is clear that this approach will not give the right exponent since to get the 
right exponent we have to choose some complicated A. Thus we shall not insist 
on getting the one-particle properties right. 



7 Electron-Phonon System 

Consider the lagrangian for phonons. 

L phonons = £ -±- ^ «,(*)«-,(*)) 04) 

Here it q = fe q + fe^ q is the phonon displacement and the phonon dispersion is 
chosen to be acoustic : O q = v s \q\. From Mahan [8] we may write down the 
hamiltonian for the electron-phonon interaction neglecting umklapp processes, 
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We choose — V e i(q) — Z\H (47re 2 )/q 2 . The overall action may be written as, 

*/3^z„ C(qn) P<\(j,nP— q<7,— n 

qer,n q<7,n qn 

+ E ^" 4 n"^ "l."""*-" + "T^ E M q/<W"q,« (97) 

qn ^ * q(T,n 

Here 77 = N ion /N e = l/Z ion is the ratio of the number of ions to the number 
of electrons and h — 2m e = 1. Also ~ 1 is the effective charge seen by 

the outer electrons. In the case of phonons the momentum transfer |q| < A/j, 
the Debye cutoff. We may integrate out the phonons first and write down the 
effective action for the electrons. 

Pqer,nP-qer,-n 

q(T,n qer,n qn 

2 P °M|n q 

qcr,f7 ,n 

In order to test for superconductivity we have to compute these quantities. The 
first is the one-particle dynamical density of states which is related to the full 
dynamical propagator. The second is the momentum distribution which should 
not have a discontinuity at the Fermi surface. The third is Yang's off-diagonal 
long range order correlation function. Unfortunately to get these right we have 
to get the one-particle propagator right. We shall now postulate that current 
algebra is more sacred than Fermi statistics and set (following the suggestion of 
A.H. Castro- Neto), 

^ slow {^a,t) ^ e^ n ^^ ^ (99) 
J(xcr, t) w -p° Vil(xcr, t) (100) 
In three dimensions we may write, 

1 ^ e 9 ^ - 1 

(p^,np-da,-n) = — T, E + q/ 2 )(! - M k ~ <l/ 2 )) TkTq" 

P k 1 m Zn 

k% V [°° J f +1 J . / 1 
-r-r I ay j a cosu 



-i/3 (2tt) 2 " J_! \ e P{vFV+^ cosO] + 1 

I \ e l3v F q cos0 _ 1 



cos8] _|_ ij iv F q cosO - z 

7T~^2 / dx Tanh 
-i/3 J (2tt) 2 J 



p^L x — (101) 

IVFq X + Z n 
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8 Conclusions 



It should be clear to the reader by now that bosonization is a powerful method 
that could hold the key to understanding a wide variety of phenomena using 
an extremely general and economical set of rules. The sea-boson method has 
reduced what was till recently, difficult research problems into difficult home- 
work problems. Now the hope is that the hydrodynamic approach will reduce 
everything to homework problems. The preferred method favored by condensed 
matter theorists involves starting from experiments and working backwards and 
trying to come up with some minimal description may sound as the most eco- 
nomical way of proceeding. However this approach is fraught with ambiguity 
and lacks the generality that the present theory has. In fact this approach of 
starting from experiments may be compared to this humorous situation. Imag- 
ine an athlete when asked to run a marathon stands on the finish line. When 
the umpire shouts 'ready!', he takes a few steps back ...'get set'. ..'go!'. The 
athlete runs like the wind and declares himself a winner of the marathon after a 
fraction of a second. Such an athlete would be the laughing stock of the sport- 
ing world. Yet it is this approach that most condensed matter theorists favor. 
The reader may argue that there is no telling where the starting line is, it could 
be 22 light years away (as in String Theory or M-theory) rather than 22 miles 
away. While this is a valid objection to my stand in principle, in practice it 
is not. Rutherford and others have already showed us using the starting-from- 
experiments approach that (for most purposes) the atom is the building block of 
matter and the atom is made of positive and negatively charged point particles 
making the atom neutral and these obey quantum mechanics. Thus there is no 
need to repeatedly ask the experimentalist for data. All that these experiments 
in condensed matter will tell us is that matter is made of positive and nega- 
tive charges and quantum mechanics is important. We already know that. We 
have to now work out the consequences of this knowledge. The main aim of 
this preprint is to convince the reader of the urgency of somewhow making this 
scheme work for then we will be able to compute the observable properties of a 
theory of neutral matter where the only adjustable parameters are the relative 
concentrations of the various elements that make up the material. 



9 Appendix A 

Here we find the most general solution to the recursion relation for the phase 
functional that determines the statistics. We find that it has to be in the specific 
form that has been used in the main text. Consider Eq.( 6) with m being in 
general, a function of the density and also a function of the pair of points x, x 
in such a way that m([p];x,x ) = — m([p];x ,x). We may formally Fourier 
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transform the phase functional as follows. 

$([p];x)= J D[P] e iP -'<X[P];x) (102) 
Therefore, Eq.( 6) may be rewritten as follows. 

0([P];x) (e- p ( x ') - 1) - 0([P];x ) ( e -* p ( x > - 1) = tt m([P];x,x') (103) 
This means that m has to be of the form, 

([P];x,x ) = (e-* p ( x ') - 1) ( e - p ( x ) - 1) (p([P];x) -P([P];x )) (104) 



m 

In other words, 

0([P];x) 0([P];x') 



e -»P(x)_ 1 e -iP(x')_ 1 

This means, 

<X[P];x) 



= P([P];x)-P([P];x) (105) 
P'([P];x) (106) 



e -iP(x) _ X 

Or, 

<&([p];x) = A([{p(y) - <5 d (y - x)}];x) - A([p];x) (107) 

Unfortunately for current algebra to be strictly obeyed, we must have A([p]; x) = 
A([p]) independent of x. This violates the statitstics requirement. Thus we shall 
have to make do with recovering the current algebra as a limit. 
Next we prove Eq.( 36). To do this we write, 

v r (x) = J D[P] e tP - p v r (x;P) (108) 

We plug this into the equation below 



to obtain, 

5 



-id-. 



J D[P] e«> v^';P)] = [-iO.j^ry j D[P] 5 -(*;P)] 

(110) 



J D[P] e lP - p \d x P(x) v r (x';P)-d x <P(x) v r (x;P)] =0 (111) 

Or, 

d x P{x) i r {x-P)-d x jP{x) v r (x;P) = (112) 

In other words, 

v r (x;P)=G([P])d x P(x) (113) 

and Eq.( 36) follows. 
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10 Appendix B 

Here we would like to ascertain whether or not the velocity operator of fcrmions 
is irrotational. To this end we note that if, 

J = -pVII (114) 

then in three dimensions, 

pVxJ = VpxJ (115) 

Thus we would like to verify whether or not this is obeyed. In general we may 
write, 

J ( x ) = -^E k 4 +q/2 ^- q / 2 ^ qx (us) 

kq 
kq 

From this we can see that, 

pVxJ-VpxJ = -d^^ [(q- q ) X k] ^ +q72 v_ q72 4 +q/2 c k _ q/2 e -(<W).x 

kq k'q' 

(118) 

This operator is not identically zero but is zero in suitably restricted Hilbert 
space. We have assumed in the text that this is the case even though it is 
strictly speaking not right. Besides, current algebra is obeyed only if the velocity 
operator is strictly irrotational as already shown. 
Consider the hamiltonian, 

q#0 q#0 

where [-Xq,Pq] = i and all other commutators involving any two of these is 
zero. There seems to be only one inequivalent way in which we may diagonalisc 
Eq.( 119). It is to decompose Eq.( 119) into oscillators with just one momentum 
label 6 q . This means we may write, 

^ = ^(b-,-bl) (120) 



p q = V^(6 q + 6 + _ q ) (121) 



q/0 q#0 

= E e « b % (l 22 ) 

q#0 
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This decomposition describes bosons. To describe fermions we need a dispersion 
of the kind vf\q\ rather than e q = q 2 /(2m). This is impossible to achieve with 
the hamiltonian in Eq.( 119). This fact can be vcrficd independently using 
the equation of motion approach that also suggests that the dispersion of the 
modes is e q . Thus we need an additional phase functional A to give us a linear 
dispersion. 
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